Transient Dynamics of Pinning 
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We study the evolution of an elastic string into the pinned state at driving forces slightly below 
the depinning threshold force F c . We quantify the temporal evolution of the string by an activity 
function A(t) representing the fraction of active nodes at time t and find three distinct dynamic 
regimes. There is an initial stage of fast decay of the activity; in the second, intermediate, regime, 
an exponential decay of activity is observed; and, eventually, the fast collapse of the string towards 
its final pinned state results in an decay in the activity with A r ~ (t p — iy , where t p is the pinning 
time in the finite system involved. 
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Driven dynamics of disordered elastic media, which has 
become a paradigm for a diversity of physical systems, 
is characterized by the existence of a depinning force 
F c , which separates the pinned state at "sub-threshold" 
forces F < F c and the sliding regime at "super-threshold" 
forces F > F c . At finite temperatures the system can 
move at sub-threshold forces as a result of thermally ac- 
tivated jumps between the disorder-induced meta-stable 
states. The corresponding activation barriers diverge at 
small forces, U(F) — > oo as F — > 0, giving rise to glassy 
dynamics with a highly nonlinear response to infinitesi- 
mal drives (creep), aging, memory, and hysteretic 
effects 0-0]. On the qualitative phenomenological level, 
the observed memory and hysteretic behaviors can be 
viewed as a result of slow relaxation of the relevant ac- 
tivation barriers upon altering the drive |^|, yet the de- 
tailed quantitative description of transient glassy dynam- 
ics remains an appealing task. The static sub-threshold 
characteristics of the string were studied in Q and criti- 
cal depinning and transient behavior in the sliding regime 
in [ [To|| . In this letter we investigate numerically and an- 
alytically the transient sub-threshold zero-temperature 
dynamics of an clastic string, taken as an exemplary sys- 
tem for a driven elastic disordered medium. We study 
the time evolution of the string as it slows and eventu- 
ally stops upon applying a constant driving force in the 
vicinity of the depinning threshold, F < F c , where the 
transient period and the distance traveled by the string 
are large and diverge as F — > F c . 

The graphic representation of a typical time evolution 
of the string velocity v(x, t) is given in Fig. [J. It demon- 
strates the avalanche-like sub-threshold dynamics of the 
string. The avalanches advance via the propagation of 
kinks along the string. We quantify the temporal evolu- 
tion by the activity function A{t) defined as the fraction 
of the moving nodes at time t, which is proportional to 
the velocity of the center of mass of the string |nj , and 
find three distinct regimes of the temporal evolution of 



the string: (i) an initial regime of a fast decay of activity, 
according to a power law; (ii) the exponential decay of 
activity, A ~ exp [— (t/t&)~ 1 ] and (iii) the final avalanche- 
like immobilization stage, which appears only for strings 
of finite length for which the activity vanishes and the 
string gets pinned at a finite time t p . Shortly before t p 
the activity function drops rapidly to zero following the 
power law, A r ~ (t p — t)^, with tp w 0.63. 

Model - We consider an elastic string at zero temper- 
ature placed on a disordered plane and subject to a con- 
stant drive F not exceeding the critical depinning force 
F c . The overdamped equation of motion is rjdy/dt = 
—STL/Sy, where y(x,t) is its lateral displacement at time 
t, rj is a friction coefficient. The string has a length L 
and is aligned along the x direction. Its energy is 



H = 



dx 



dy 
dx 



o h£ +U(y(x),x) 



with the string elastic constant C . The string and the 
media are periodic in x with periodicity L. 

In the simulations we use a discrete version of the con- 
tinuum equations. The x direction is discretized by a set 
of lines (rails) of equal spacing, with unit length, con- 
straining the motion of discrete string elements (nodes). 
Energies are measured in the units of the elastic con- 
stant C . The pinning sites are randomly distributed on 
the rails according to a specified density p — 0.1. The 
pinning sites are triangular wells with unit width and the 
depth/strength U p = 0.1. Forward Euler time marching 
is used to advance the system in time. The initial state 
for the string is taken as a straight line. The statistics 
were collected on 30 samples of length L — 8192. 

Results - To link to past research, we measured the 
geometric characteristics of the subthreshold behavior. 
The roughness of the string at time t is Wt(x) = ([y(x + 
x',t) — y(x' , t)] 2 ) 1 / 2 , where (. . .) denotes a positional av- 
erage over x' and a disorder average. The roughness in- 
creases with x (for x <C L) until it saturates at x ~ £m (t) 
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with w t {x) » £±(t) (see Fig. ||). Correlation lengths 
increase with time and reach the final values £|| and £j_ 
in the pinned state (t > t p ), where the roughness fol- 
lows a power law ~ x^ for 1 <C x <C £|| Qj. We 
find the roughening exponent £ to grow linearly as the 
applied force approaches its threshold value F c (see in- 
set of Fig. ||). This behavior allows for an alternative 
determination of F c resting on the analytical estimate 
( = 1 at the transition p3[ . This contrasts with the 
usual procedure where two parameters F c and £ have 
to be simultaneously determined as fitting parameters. 
The criterion ((F c ) = 1 yields as an estimate for the 
critical depinning force F c — 0.3016, which is in excel- 
lent agreement with our value of the critical force deter- 
mined independently from the double-logarithmic plot of 
roughness (F c = 0.3105) and also with the earlier finding 
Q(F C = 0.3058). Note a discrepancy with the results in 
Rcf. [(To), where £ = 1.25 at the depinning transition has 
been found in conflict with the analytical result of [fl3|| . 
The reason for this discrepancy is not yet clear. 

Next, we studied the behavior of the pinning distance 
D p , which is defined as the ensemble averaged distance 
that the center of mass of the string traveled before be- 
ing completely pinned D p = (1/L) Yh=i \Vi(tp) - 3/i(0)|. 
We verified the scaling behavior of D p ~ (F — F c )~ a 
(cf. Fig. |). For the given string length L = 8192 the 
finite-size effects become negligible for the force depen- 
dence. We found a = 2.56 ± 0.08 for F c = 0.3016 in 
good agreement with the results of ^ for the ID CDW. 
The other relevant quantities £|| and determined from 
the saturation of the roughness also scale with exponents 
grouping around a w 2.4. The results are summarized 
in the Table Q, which also gives the values of the critical 
force determined independently from the divergences of 
the quantities involved. 

Having established the geometric characteristics of 
subthreshold dynamics, we now turn to a detailed study 
of the temporal evolution of transient behavior. We char- 
acterize the temporal evolution of a finite system toward 
the pinned state by the activity function A{t) defined as 
the fraction of the moving string nodes moving at the 
given instant t. For a discrete system evolving in time 
by steps 5t, we define a node as active at time t if that 
node has moved at least one lattice unit in the previous 
100 time steps. The activity function is then the ensem- 
ble average of the fraction of active nodes. In Fig. ^ we 
illustrate the temporal behavior of the activity function 
for the representative subthreshold force F — 0.29. 

It is important to quantify not only the fraction of ac- 
tive nodes, but also the spatial structure of pinned and 
moving nodes, which can be quantified by the charac- 
teristic length ipi n of pinned segments (connected re- 
gion of pinned nodes). Initially, L v - ln will be of the or- 
der of unity, since the string position is not correlated 
with the location of the pinning centers and in the final 
pinned state L p ; n = L (see inset of Fig. 0). The initial 



transient regime (i) is expected to last until t% given by 

£pin(*«)=£|| II- 

At t > relatively large segments of the string are 
pinned, and the remaining activity is due to the motion 
of mobile segments over rare regions, where the pinning 
centers are underrepresented. These mobile segments will 
advance until the driving force is balanced by other pin- 
ning centers or by the string tension. We now construct a 
lower estimate for the velocity of an infinitely long string 
as a function of time. We expect the finite string in 
regime (ii) to follow this dynamics, before finite-size ef- 
fects set in in regime (iii). We divide the rails in ele- 
ments of unit length. Then a finite fraction p of elements 
is free of a pinning force. We overestimate the effect of 
pinning by assuming that the other elements contain an 
infinitely strong pinning force. Let us now consider a 
string segment that is pinned on two rails of distance 
say at y(0) = y(l) = 0. If the rails in between are free 
of pinning centers, the segment relaxes into the configu- 
ration where the driving force is balanced by the string 
tension. A simple estimate shows that the front of the 
segment will then reach a distance yi ~ Fl 2 /C after a 
time ti ~ yiij/F ~ nl 2 /C. Before this time the ends of 
length x(t) ~ y/ Ct/rj of the string segment come to rest, 
and the center of mass of the segment moves with a ve- 
locity vi (t) ~ {F/rj) [1 — yj Ct/rjl 2 ] . During this relaxation 
the string segment moves over an area Ai ~ Fl 3 /C. The 
relaxation takes place in this undisturbed way only if all 
the rail segments in this area are free of pinning, which 
occurs with probability pi ~ p Al . 

To find a lower bound to the string velocity we di- 
vide the whole string into segments of arbitrary length 
I, which serves as variational parameter. Although the 
string can be taken as unpinned initially, the velocity is 
underestimated by formally considering the links between 
these sections as pinned. A fraction pi of such segments 
relaxes without encountering pinning rail elements and 
contributes ~ pivi(t) to the average string velocity. This 
contribution is a lower bound to the velocity for all pos- 
sible I > yj Ct/rj, and the best estimate is obtained by 
v min (t) ~ maxipivi(t) ~ p^/ c )( ct / ? ') 3/2 , which is finite at 
all times. Consequently, the activity function is also ex- 
pected to decay not faster than A(t) > exp(— const. i 3 / 2 ). 

In the finite system this exponential decay stops when 
the characteristic distance between the moving segments 
Lpi n ~ / /pi becomes of the order of the system size and 
the string experiences a deficiency of the moving nodes. 
Our conjecture is that the origin of this "superfast" stop- 
ping regime is the avalanche-like disappearance of the 
moving nodes of the string. Once pinned, a node cannot 
depin any more; moreover, it causes the immediate pin- 
ning of the neighboring nodes, in a close analogy with the 
avalanche clustering at the final stage of the coagulation 
process |l5[ . This analogy suggests the power-law termi- 
nation of the motion according to A r (t p — t) ~ \t p — t\^, 
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where t p is the time of the total immobilization. Al- 
though t p is subject to strong sample-to-sample fluctua- 
tions because it is determined by fluctuations on the scale 
of the system size, the final regime can be recognized only 
by averaging the activity considered as a function of the 
reverse time t p — t; that is, the final regime can be seen 
only if the averaging over N samples is performed accord- 
ing to A r (t p —t) = (1/N) J2iLi Ai(tp — t) in contrast to 
the usual average A(i) = {1/N) Y^ =1 A^t). 

We give lower and upper limits on the exponent if). 
Consider the process in which the node that came to rest 
enhances the probability that the nearest node will stop 
at the next moment of time. In this case the boundary 
between localized and mobile nodes will move some aver- 
age velocity, and the activity will decrease 
which means that the upper limit on ip is V'uppcr = 1- 
This type of behavior can be expected when F is far be- 
low F c . Another limiting type of terminal dynamics oc- 
curs if immobile nodes do not affect dynamics in mobile 
regions. In this case the boundary between the immo- 
bile and mobile regions moves randomly, and the finite 
string has a nonzero chance of stopping when the bound- 
aries enclosing the moving region collide with each other. 
Considering this process in reverse time, we conclude that 
A r oc (tp — t) 1 / 2 , and we get a lower limit Slower = 1/2- 
Note, that the possibility of fragmentation of the moving 
segment of the string into several subsegments does not 
affect this estimate. In the more general case of a moving 
manifold of internal dimension d (d = 1 for the string) , 
similar arguments result in d/2 < ip < d. 

Numerical simulations clearly show distinct relax- 
ational regions. After the initial relaxation (regime (i)), 
the activity crosses over to the exponential decay regime, 
A(t) ~ exp [— (i/iA) 7 ], 7 — 1. The precision of our 
data is not sufficient to determine a more precise esti- 
mate of 7; from the above consideration we could expect 
1 < 7 < 3/2. The relaxation time t& diverges as the ap- 
plied force approaches the threshold value (see Fig. ||). 
However, although one could expect that t& — t% and 
therefore that t A ~ \F C -F\"* (since £ ~ \F — F C \"), with 
zv = 4/3 using the exponents of the depinning transition 
fiof , we have found a different exponent: t A ~ \F C — F\ A 
with A = 2.74. A careful examination of a final stage of 
the relaxation process reveals a very sharp drop of the 
activity function just before the complete immobiliza- 
tion of the string. To analyze this stage, we introduce 
the reverse time representation: we present the activ- 
ity function as a function of t r = t p — t, counting time 
backwards from the moment of achieving final immobi- 
lization. Using this representation enables us to carry 
out an ensemble averaging despite the fact that different 
samples have different pinning times. The plot of hiA r 
vs. ln(t p — t) reveals a power-law final behavior (see Fig. 
g): Ar~(tp- t)^, with ip = 0.63 ± 0.06, which is within 
the above upper and lower estimates. We note that the 



reverse time plot of In A r reconfirms the existence of the 
intermediate regime with the essentially same character- 
istic time tA as was found from the direct time picture. 



In summary, we investigated subthreshold dynamics of 
a pinned elastic string and identified three distinct tran- 
sient regimes: (i) a fast initial relaxation, (ii) an inter- 
mediate exponential decay of the activity resulting from 
residual motion in the exponentially rare regions free of 
defects and (iii) a novel avalanche-like terminal relax- 
ation to the pinned state, A r ~ (t p — t)^ , resulting from 
the finite-sized effects. This final stage exhibits striking 
similarity to coagulation dynamics. 
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FIG. 1. History of the string velocity (gray-scale plot of 
v(x, t) using white for zero velocity and black for large veloc- 
ity). Note that in the active regions (avalanches) the velocity 
of the beads reaches the same values at early and late stages 
of the evolution, which differ only in the number and size of 
avalanches, which is measured by the activity. 



FIG. 4. Time dependence of the activity A for the force 
value F = 0.29 in linear and semilogarithmic representation. 




FIG. 2. Plot of the roughness w t {x) at intermediate time 
(t < t p , truncated for large x beyond saturation) at final stage 
(t > t p ), where it saturates on the scale of the correlation 
lengths £||,±. The nonmonotonic shape of Wt(x) is an artifact 
of the finite number of samples used for disorder averaging. 
Inset: time evolution of the typical size L p i n of pinned string 
segments. 
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FIG. 5. Plot of the activity A r as a function of the reverse 
time (t p — t). The dashed lines represent linear fits to the 
regimes (ii) and (iii) introduced in the text. 
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FIG. 3. Inset: The roughening exponent £ as a function 
of the driving force F (stars: measured values, dashed line: 
linear fit), which assumes the value £ = 1 at the critical force 
F = F c . Main plot: force dependence of the pinning distance 
Dp (crosses) and of the relaxation time 4a (triangles) with 4 
linear fits (dashed lines) in the double-logarithmic represen- 
tation. 



TABLE I. Values of the scaling exponent a and the critical 
force F c as determined from the divergences of D p , £±, and 
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